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ABSTRACT 


We color the vertices of a graph G, so that no two adjacent vertices have the 
same color. We would like to do this as cheaply as possible. An efficient coloring would 
be very helpful in optimization models, with applications to bin packing, examination 
timetable construction, and resource allocations, among others. Graph coloring with 
the minimum number of colors is in general an NP-complete problem. However, there 
are several classes of graphs for which coloring is a polynomial-time problem. One 
such class is the chordal graphs. This thesis deals with an experimental algorithm 
to approximate the chromatic number of an input graph G. We first find a maximal 
edge-induced chordal subgraph H of G. We then use a completion procedure to add 
edges to H, so that the chordality is maintained, until the missing edges from G are 
restored to create a chordal supergraph S. The supergraph S can then be colored 
using the greedy approach in polynomial time. The graph G now inherits the coloring 


of the supergraph S. 
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DISCLAIMER 


The computer programs in the Appendices are supplied on an “as is” basis, 
with no warrantees of any kind. The author bears no responsibility for any conse- 


quences of using these programs. 
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Ai. INTRODUCTION 


Large-scale scheduling and timetable problems arise in many activities, from 
crew scheduling on airlines to scheduling of classroom periods and teachers at Uni- 
versities to register allocations in a computer CPU. These are often modeled as graph 
coloring problems that are then subjected to a variety of strategies for solution by 
computer. For example, consider the storage problem at military ammunition supply 
points, where certain types of ammunition cannot be stored in the same bunkers. 
What is the minimum number of bunkers needed at an ammunition supply point to 
ensure the safe storage of all types of ammunition? We use a graph-theoretic model to 
formulate this problem. We construct a graph G, with vertices representing ammuni- 
tion types and where the existence of edges between vertices represents their pairwise 
incompatibilities. This problem reduces to that of finding an optimal coloring of the 
vertices of G. These types of optimization models have applications in bin packing, 
examination timetable construction, and resource allocations. Unfortunately there 
is no known algorithm for this problem which will predictably provide an optimal 


solution in a reasonable time. 


A. DEFINITIONS 

The graph terminology used in this thesis can be found in most textbooks on 
graph theory. For undefined terms and notation see Bondy and Murty [Ref. 1], or 
West [Ref. 2]. 


1. Undirected Graph 

A simple undirected graph G = (V, E) consists of a vertex set V(G) = {u,..., Un} 
and an edge set E(G) = {e1,...,;€m}, where each edge is an unordered pair of ver- 
tices (see part (a), Figure 1). We use uv to denote the edge {u,v}. When we have 
uv € E(G), then we say that “u is adjacent to v” and “v is adjacent to u”. The 


vertices of an edge e are its endpoints. If the endpoints are the same this edge is 


e a =) a e ra 


G=(V,E) G=(V,E) Gy=(V,E") 
n=5, M=/ E’ =(ab,bc,ce,de) A=(b,c,d,e) 


Figure 1. (a) Undirected graph G of order n and size m, (b) A subgraph H with edge 
set £7’, (c) An induced subgraph with vertex set A. 


considered a loop. In a graph G, a list of vertices [vo, v1, v2,..., v1] is a path of length / 
from vertex vp to wy if {v;-1,u;} € E(G), for all2 =1,2,...,/. A path is called closed 
if vp = v;. A cycle is a simple closed path of distinct vertices. A graph G is connected 
if between any two vertices there exists a path in G joining them. The nezghborhood 
of v is the set of adjacent vertices Adj(v) and is denoted N(v). We consider only 


connected, loopless, finite graphs without multiple (redundant) edges. 


2. Miscellaneous Graphs 

The order of a graph G, denoted |V]|, is the number of vertices in G. We use 
|E’| to denote the number of edges, or the size of G. The complete graph of order n is 
denoted A, (see part (c), Figure 2), and is an undirected graph in which every pair 
of vertices is adjacent. The complement of a graph G is the graph G = (V, F), where 
E={(z,y)€ Viz # y and (z,y) ¢ E} (see part (b), Figure 2). A subgraph of Gis a 
graph H, such that E(H) C E(G), V(H) C V (G)(see part (b), Figure 1). If ACV, 
the A-induced subgraph of G is the graph Gy, such that E(G4) = {ry € E(G) |x € 
V (G4) and y € V (Ga)} (see part (c), Figure 1). 


3. Graph Parameters 
For the following definitions let G = (V, £) be an undirected graph. A coloring 
of a graph is an assignment of “colors” to the set of vertices V so that adjacent 


vertices have different colors. The symbol y(G) denotes the smallest number of colors 


to 


e a e a e a 


G=(V,E) G=(V,E) Ks 


Figure 2. (a) Undirected graph G, (b) Complement graph G with edge set EF, (c) 
Complete graph Ks. 


required and is called the chromatic number of the graph G( see figure4). A subset 
V (H) © V(G), with |H| = q, is a q-clique if it induces a complete subgraph Kj. 
A clique H is mazimal if there is no clique of G which properly contains H as a 
subset and H is maximum if no other clique in G has a larger order. The number of 
vertices in a maximum clique of G is called the clique number of G, and is denoted 
w(G). An independent set S C V (G) is a set of vertices in which no two vertices are 
adjacent. The number of vertices in a independent set of maximum order is called 
the independence number of G, denoted a(G). A subset K © V(G) such that if 
ry € E(G), then x € K or y € K, is called a verter cover. A vertex cover for G is a 
set of vertices that is collectively incident to all the edges in E(G). The vertex cover 


number of G is the number of vertices in the minimum vertex cover, denoted @ (CG). 


4. Chordal Graphs 
A graph G is chordal (or triangulated) if it does not contain any cycle of length 


greater than three as an induced subgraph(Ref. 3]. An ordering of the vertex set V, 
with n = |V|, is {vo, v1, v2,..-, Un}. A vertex v is simplicial in G if N (v) is a complete 
subgraph. If 6 is an ordering of V such that each v; 1s a simplicial vertex of the induced 


v,, then it is called a perfect elimination ordering (see Figure 4). A 


subgraph G,, 


peers 


successor of v; with respect to the ordering 6 is a vertex v; € N(v;), where 2 < J, 


and is denoted Suc(v;)[Ref. 4]. A chordal graph G is of the class of perfect graphs, 


in Which w(G4) = y (Ga) for all AC V(G). The class of perfect graphs have been 
extensively researched since the 1960’s by well-known mathematicians. Golumbic’s 
book “Algorithmic Graph Theory and Perfect Graphs” [Ref. 3] does an excellent 
job in explaining the difficult practical problems related to the structure of perfect 


graphs. 


5. NP-Complete Problems 

The class of problems with complexity bounded by a polynomial in the size 
of the input is denoted P. We consider a problem solvable in polynomial time by 
a deterministic algorithm as being tractable. That 1s, for each input of size n the 
worst-case running time is O(n*) for some constant k. We define problems that 
require superpolynomial time as being intractable. 

The class VP contains those decision problems that are “solvable” by a non- 
deterministic polynomial-time algorithm. Such an algorithm, tn a sense, tries all 
possibilities simultaneously, applying polynomial-time computation to each guess in 
parallel. This type of algorithm should not be confused with a parallel implementation 
of a deterministic algorithm. If any of the computations results in a yes or possibly a 
no, then the algorithm is a success. The algorithm is successful if it works, even if the 
answer to the current decision problem is negative. The non-determinism concerns 
the multiplicity of paths, and not whether the search is successful. It 1s easy to see 
that if we can do many computation paths in parallel and one of these is completed 
in polynomial time, then we can do that one alone in polynomial time. Therefore 
Bac P. 

Most mathematicians believe that the classes P and NP are different classes, 
although it has not yet been proven that P # NP. The class P, loosely, consists of 
those problems that can be solved quickly, while the class NP consists of problems 
for which a solution can be verified quickly. We define a problem X as NP-hard if 
every problem instance in NP can be reduced to an instance of X in polynomial 


time. A problem is NP-complete if it is in NP and is NP hard. We now have the 


class of P in the class of NP, but the problems that are NP-complete are also in 
NP, soa problem in P is probably not N P-complete. Since no one has come up with 
a polynomial algorithm for a NP-complete problem, thus proving that P = NP, 
we can assume the intractability of NP-complete problems. Problems that on the 
surface seem no harder than sorting, graph searching, or network flow are in fact 
NP-complete. Thus, it is important to become familiar with this class of problems. 


For a detailed discussion of N P-completeness see Cormen, Leiserson, and Rivest[Ref. 


5] or Garey and Johnson|Ref. 6]. 





If. EXISTING ALGORITHMS 


Ore’s book “The Four-color Problem” [Ref. 7] shows that considerable litera- 
ture in the field of graph theory deals with the coloring of graphs. Many algorithms 
exist for graph coloring, but there are only a few fundamentally different approaches 
to the problem. Finding the exact chromatic number of a graph is a NP-complete 
problem, so it 1s no surprise that no fast algorithm exists. We will explain two such 


algorithms and it should then be obvious to the reader why they are impractical. 


A. EXACT ALGORITHMS 
1. Brute-Force Coloring 


The different ways to color a graph are not unique. There may be several 
proper colorings of graph G using x(G) colors. There also exist many improper 
colorings of that same graph. How many different colorings of G are there? A 
coloring of the n vertices in V using a palette P of order k is a mapping f: V — P. 


If we use all of P, the mapping is onto. Then there are 


()e-(h)ecire(etaacares een 


different colorings, and if k > y(G) one or more of these colorings is proper. We could 
easily program a computer to produce these colorings and for each, check to see if 
the coloring is proper. But the cost of computation would be exorbitant. Using the 
formula above we see that the number of colorings gets outrageous very quickly. For 
example, let n = 10, then if k = 2 there are 1022 different colorings. If k = 3 there are 
55,977 colorings and if k = 4 there are 818,521 different colorings. If we wanted to 
color a large graph in this way, say n = 100, and k = 2 there are 1.267 * 10°°different 


colorings. 


2. Greedy-Backtracking Coloring 

Given G, this algorithm finds \(G) if given a starting palette of fewer than 
\(G@) colors. It performs an extensive routine of trying to color the graph with the 
palette given, and if it determines that the palette is not large enough then another 
color is added to the palette. Eventually enough colors appear on the palette to 
successfully color the graph. 

The most significant problem with this algorithm is to determine the starting 
number of colors k. If & is too low, the algorithm consumes too much time backtrack- 
ing and recoloring the graph and if k is greater than y then algorithm will greedily 
color the graph with the available colors & and the resulting coloring may not be op- 
timal. The clique number w(G) is clearly a lower bound on y(G), since the vertices 
of the largest complete subgraph of G must all have separate colors. In small graphs, 
w(G) is very close to x(G), but in larger graphs the difference can grow significantly, 


as 1s shown in the next section. 


B. APPROXIMATION ALGORITHMS 


Several algorithms exist to approximate the chromatic number of an arbitrary 
graph. Some are better than others. Most of these algorithms deal with finding the 
largest independent set of vertices or the vertices of a maximal clique. To find the 
largest independent set is a NP-complete problem, but to find a large one can be done 


in polynomial time. 


1. Independent Sets 

The approximation algorithm for the independent set problem is based on two 
assumptions: (1) a vertex of high degree is harder to color than a vertex of low degree; 
(2) coloring many vertices with the same color is good. Recall that a(G) = w(G) is 
the order of a largest independent set of the vertices in G. Clearly the vertices in 
an independent set can be colored with the same color, therefore y(G) < n/a(G), 


where n = |V|. In small graphs, n/a(G) tends to be smaller than w(G), but in larger 


graphs n/a(G) is a much closer lower bound to ,(G). 

Let f(n, p) = E(q@) bea function to find the expected value a(G) for a random 
graph G of order n and edge probability p. With this information we can estimate y. 
In principle we can find an independent set of order E (a), and delete it. Now we have 
a graph of order n — f(n,p) = n,. We continue on the same way until E (@~) = 0 and 
thus y(n, p) =2+xX(n:—f (ni, p),p) = n/a(G), for: = 0,1,2,.... Note: o(G) = a(G) 
for graphs with edge probability p = .5 and in particular x(1000,.5) = 85, but 
a(1000,.5) = 15 and n/a(i000,.5) = 67. This points out a very peculiar problem 
with graphs of large order and density. That is to say, w(G) is a very poor estimator 


while n/a(G) is a good estimator of y(G) for graphs of large order and density. 


2. Vertex Covers 

In a graph G, if a set A C V(G) is a vertex cover then there are no edges 
in A, an independent set. Thus any minimum vertex cover is the complement of 
a maximum independent set, and so a(G) + 8(G) = n, where n is the order of 
G. The vertices in a minimum vertex cover are the only vertices considered in the 
coloring problem; all other vertices are part of the independent set and require only 
one additional color. Therefore G(G) +1 > x(G). Another upper bound on x (G) 
that is worth considering is A(G) +1. Both errors can be relatively large, though. 

The vertex cover problem is known to be N P-complete (see Cormen[Ref. 5]). 
Nevertheless, there exist good algorithms to find a vertex cover that is near optimal. 
One such algorithm uses a set C, initially empty, and a set F containing the edges of 
a graph G. We pick an arbitrary edge ry € E and add the vertices x and y to C, then 
delete any edges in EF covered by x or y. We pick another edge in / and continue 
this procedure until & is empty. The computational complexity of this algorithm is 
O(m), where m is the size of G. The vertex cover produced by this algorithm is C, 
which is at most twice the size of the optimal cover C*. Let A be the set of arbitrary 
edges picked in the algorithm. No two edges in A share an endpoint, since all incident 


edges to the endpoints are deleted before the next edge is picked. Therefore when 


two vertices are added to C, |C| = 2|A|. Any vertex cover of A must contain at 
least one endpoint of each edge in A. Since no two edges in A share an endpoint, no 
vertex in the cover is incident on more than one edge in A. Therefore, |A| < |C%|, 


and |C| = 2|Al < 2|C*|. 


3. Maximum Clique 

This algorithm, by Balas and Yu [Ref. 8], is a chordal subgraph approach for 
finding the maximum clique problem. It has two main subroutines. The first algo- 
rithm generates a maximal triangulated induced subgraph H of an arbitrary graph G 
In a computational complexity of O(n +m). The second finds the minimum color- 
ing of H, using the cardinality & of the maximum clique; then appends vertices to H 
while maintaining its chromatic number, until the resulting graph becomes a maximal 
k-chromatic induced subgraph F' of G. If F = G, we are done, since the maximum 
clique in H was also maximum in G. Otherwise we branch to subproblems consider- 
ing any clique larger than the current one must contain one of the vertices in V(G) 
but not in V (F’). We now apply the same procedure above on the new subproblems, 
each defined on a vertex set contained in the neighbor set of v € V(G)/V (F). For 
the results of this algorithm and the different variations applied to it see Balas and 


Yul[Ref. 8]. 


4. Minimal Weighted Coloring of Chordal Subgraphs 

The algorithm in the previous section was modified by Balas and Xue [Ref. 9] 
to find the minimum weighted clique and thus a coloring of a chordal subgraph H of 
G. They then extended the algorithm to include an ordering 7 of the vertices V\H. 
This ordering 7 is used to add remaining vertices to the correct color class until a 
maximal induced subgraph F' with the same minimum weighted coloring of H results. 
The final step is to modify the branching rules described by Balas and Yu [Ref. 8] to 
include the minimum weighted coloring and define the subproblems to reapply in the 


algorithm until it finds the maximum weighted clique. 


LO 


5. Edge-Maximal Chordal Subgraph 

The algorithm by Xue [Ref. 4] involves n iterations. Each iteration augments 
the partial perfect elimination ordering and adds a vertex, together with some edges, 
to the partial chordal subgraph. What sets this algorithm apart from algorithms 
like that of Dearing et al.[Ref. 10] are the way in which it chooses the next vertex 
to add into the partial chordal subgraph, and the way in which it chooses the first 
successor of a given vertex. We use a greedy approach in yon instances. During an 
iteration, we call the vertices in the partial perfect elimination ordering labeled and 
the rest unlabeled. Let U be the set of unlabeled vertices and H = (V(H), E(#)) 
be the partial chordal subgraph. For every unlabeled vertex v € U, we assign v 
a temporary first successor t(v) and a label s(v), where t(v) € V(H). The label 
s(v) =1+ |N(v)N Suey (t(v))| is the maximum number of edges that can be added 
into H if v is added into H next with t(v) being its first successor. We initially set 
t(v) = @, no temporary first successor and s(v) = 0, Vv € V. We choose the next 
vertex v € U to label and add to H such that s(v) = max{s(u)|u € U}. Ties go to 
the vertex v with the largest degree. We add all the edges to ¢(v) or to a vertex in 
Sucy (t(v)). We update t(u), s(u): For all ue N(u) NU, let r, be the number of 
neighbors of u in H that are either v or a successor of v. If ry < s(u), set v as the 
temporary first successor of u and update s(u), i.e., let r, = 1+ |Sucy (v)N N (u)]. 


ieee Su), set (2) =v, s(u) = 7r,. 
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III. A SUPERGRAPH HEURISTIC 


The basic idea is, given G, to find S, a minimum cardinality chordal supergraph 
of G. We could then color S, and let G inherit the result. What makes this a topic 
worthy of research is the fact that the first step is N P-complete. We get around this 
by finding a minimal chordal supergraph of G, ideally containing as few edges of G 
as possible. 

In this chapter we introduce and explain an experimental scheme for coloring 
an arbitrary graph. We have developed a basic algorithm which will be modified in an 
attempt to improve its performance. Each experiment consists of the generation of 100 
random graphs of order 7 to 20, to which the algorithm is applied. An implementation 
in the MATLAB programing language is given in the Appendices for all functions 


whose names appear here in the verbatim typeface. 


A. THE IDEAL ALGORITHM 


The ideal supergraph algorithm for the coloring an arbitrarygraph G consists 
of finding the closest edge-induced chordal supergraph S of G. We color S with a 
simple greedy coloring algorithm, capitalizing on the perfect ordering of the chordal 
graph S, and then let G inherit the coloring. At a minimum we have y (5) > x(G) and 
if only a few inconsequential edges were induced to find S then possibly y (S) = x(G). 
The following algorithm depicts the ideal coloring algorithm. 


Ideal Algorithm 
Input: Graph G 
Ouput: Optimal coloring of a minimal chordal supergraph S$ 

begin 

S := minimum edge-induced supergraph of G; 

VS =—aGOlorimg Ol 9: 

H inherits the coloring ¥; 

end; 


S | 
= | 
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Figure 3. The Ideal Algorithm starts with an arbitrary graph G, computes a maximal 
chordal subgraph H, and then builds a supergraph S while minimizing p. (@ = (7)) 


The problem with this algorithm is that of finding the chordal supergraph S. 
We propose, as depicted in Figure 3, to find a maximal chordal subgraph H of G and 


then find a chordal supergraph S minimizing the number of edges p. 


B. THE BASIC ALGORITHM 


This algorithm, implemented in the function project1.m, is the basis for the 
improved algorithms which follow. Generally speaking, it produces a random graph, 
ensures the graph is connected, and computes its chromatic number(see figure 4). It 
then finds a maximal chordal subgraph of the random graph and a minimal chordal 


supergraph of both graphs. Lastly it computes the relative difference between the 
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Figure 4. The Basic Algorithm starts with an arbitrary graph G, computes a maximal 
chordal subgraph H, and then builds a supergraph S$. Graphs H and S vertices are 
in the perfect elimination ordering 0. 


chromatic numbers of the supergraph and the original graph. 


1. Random Graph Generation 

The basic Algorithm A, step 1, generates a random input graph G = (V, E) 
of order n and size m, with vertices labeled as v,...,u,. A random sparse adjacency 
matrix is generated by the function unigraph.m and is used to represent the undi- 
rected graph G. An edge in G exists between the vertices v; and v;, represented by a 


1 in the (2,7) entry of the adjacency matrix, with probability p, where 0 < p < 1. 


2. ‘Test for Connectedness 

The function to ensure that a graph is connected uses a depth-first search 
algorithm and has computational complexity of O(n + m) as discussed in Roberts 
[Ref. 11, page 445]. For a detailed discussion of depth-first search see Tarjan [Ref. 
12]. Since a graph of order n has at most (3) edges. we have O(n7) steps. In this 
function each time we traverse down a path to the end and return to the beginning 
without visiting every vertex in the graph we identify a connected component of the 
graph. When a disconnected graph is discovered we discard the graph and return to 
step one. 

The function connect.m starts with a list v containing the initial vertex v; 
and visits each vertex using depth-first search. At each vertex v; there are three 
possibilities: (1) there are adjacent vertices which have not been added to the list of 


vertices v. in which case we pick the lowest-indexed vertex. add it to the list v, delete 


the edge and continue the search from the new vertex; (2) there are adjacent vertices, 
but all of them appear in the list v, thus we delete all the edges and backtrack on the 
list v until we find an adjacent vertex and continue the search; or (3) there are no 
adjacent vertices, in which case we backtrack on the list v as described in case two. 
If in the process of backtracking to find an adjacent vertex we end up at the initial 
vertex v, and |v| #4 n, then we have found a connected component in a disconnected 


graph. 


3. Maximal Chordal Subgraph Computation 

A maximal chordal subgraph H of the input graph G is computed in step 
two of the basic algorithm. H is found by using an algorithm of Dearing, Shier and 
Warner [Ref. 10]. This is a polynomial algorithm used in optimization problems to 
solve large systems of linear equations. The algorithm has worst-case time complexity 
O(mA), where A denotes the maximum vertex degree in G. The mchord.m function 
is given a starting vertex v,; and the sparse adjacency matrix for G, and produces 
a list of vertices denoting a perfect elimination ordering and a list of edges for the 
maximal chordal subgraph found. 

The function mchord.m generates a perfect elimination ordering list v, starting 
with the first vertex given. It maintains a list s of all the vertices of G not in v and 
a list E of all the edges in the maximal chordal subgraph. It also builds an incidence 
matrix N of order nxn containing a | in the (2,7) entry if vertex 2 is adjacent to 
vertex ) when both vertices have been considered in the maximal chordal subgraph. 
The matrix N is initialized to all zeros and list & is empty at the beginning of the 
procedure. The next step in the mchord.m function is to generate a loop to visit each 
vertex in graph G. The current vertex is denoted v9. Step one in the loop is to find 
all the vertices adjacent to vo, using a function adj-.m. 

The function adj.m is given the adjacency matrix G and a vertex vg. The 
function determines the vertices adjacent to vg by examining the ith row, representing 


vo, of the matrix G and returning the indices of all columns 7 which include a | in 
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the (7,7) position. The function returns a list u of column indices representing the 
adjacent vertices and a list FE of edges incident to vg. 

For each vertex u that is adjacent to the current vertex uo, if N(u) GC N(vpo) 
then N(u) := N(u)U {uo} and F := EU {u, vo}. In other words, if the neighborhood 
of u is a subset of the neighborhood of vg we increment the value for the vertex u in 
the set N and add the edge to &. At this time we eliminate the edges which have 
been added to the subgraph from the adjacency matrix G, so that they will not be 
considered again later. 

From the set N we choose the vertex with the largest value to become our new 
Vg and we add it to the list of v and eliminate it from the list of s. Now we repeat the 
procedure until all the vertices in G are added to the list v. The result is the reversal 


of a perfect elimination ordering of G. 


4. Missing Edges 

The functions mkadjmat.m and mkedges.m are used to derive an adjacency 
matrix from a list of edges and to create a list of edges from an adjacency matrix, 
respectively. When we subtract the adjacency matrix H from the adjacency matrix 
G, we produce an adjacency matrix F representing the edges of G missing from H. 
The list of missing edges is denoted em = ({v;, v;}|v; is adjacent to v; and v;,v; € F). 
Anytime we reorder the vertices of H we must translate the list of missing edges em 


into the new ordering using the function trans.nm. 


5. Maximum Clique 

The clique number w(H) = y(#H), since H is a chordal Perfect graph. It 
is trivial that y(H) < x(G). Therefore, we use w(H) as a lower bound on y(G), 
because it is the easiest to find. The algorithm to find the maximum clique of a 
chordal graph by Gavril [Ref. 13], has a computational complexity of O(n +m). We 
use this algorithm in the function mclique.m to determine the w(H). We maintain 


a list S which holds the number of times the lowest-indexed vertex was a member 
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of a previous clique. As we visit each vertex in the perfect elimination ordering we 
determine its neighbors and store them in z. If x = @ then the current vertex has no 
neighbors and it is its own maximal clique. If 2 4 @ then we only need to consider the 
neighbors of higher index in the perfect elimination ordering since all vertices of lower 
order would have been eliminated and not considered in finding the next maximal 
clique. We store these vertices in X. We now update the value in S' for the smallest- 
indexed vertex u in X by S(u) = max{S(u),|X|— 1}. If the number in S of the 
current vertex is less than the current maximal clique, we output the maximal clique 
X and update the maximal clique number if the current maximal clique is larger. If 
the number in S is equal to or greater than current clique number, then the maximal 
clique has already been identified and we continue. In layman’s terms, we visit each 
vertex in the perfect elimination ordering, cutting off the portion of the graph we just 
visited. We then look forward, relative to the ordering, and determine the maximal 
clique. If it is larger than the current clique number we update the clique number 


and move on until we reach the end of the perfect elimination ordering. 


6.  Greedy-Backtracking Coloring Scheme 

The greedy-backtracking coloring scheme uses an algorithm defined by Bender 
and Wilf in [Ref. 14]. They give a detailed analysis of the run-time complexity of 
this algorithm on arbitrary graphs. The idea is to visit each vertex in the order given 
and determine which of its neighbors have been colored. We always start with the 
cheapest (or lowest) color and, having greedily colored the first & vertices, find the 
cheapest available to color the current (& + 1)st vertex. This algorithm must be given 
a palette of available colors, and if in visiting a current vertex we run out of available 
colors, we simply backtrack to the last vertex colored and determine if it is possible 
to increase the color to a (the next highest). If it is possible to increase the color we 
do and continue our search. If it is not possible we must backtrack further, erasing 
the current coloring scheme until we find a vertex color which can be increased. If 


in the backtracking process we return to the initial vertex, then we did not supply a 
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sufficient number of colors in the palette. In this case we must increase the number of 
colors available on the palette and try the procedure again. The algorithm terminates 
with a coloring that uses the least number of colors necessary for an proper coloring 


scheme of the graph. This least number of colors is of course y(G). 


7. Maximum Cardinality Search 

The vertices of chordal subgraph H are ordered in the original perfect elim- 
ination ordering returned by the mchord.m function and, consequently, most of the 
missing edges are at the front of the ordering. It would be more beneficial for our 
algorithm to have the missing edges towards the back of the ordering, since these are 
the vertices considered early in the edge completion scheme. We use the maximum 
cardinality search (MCS) algorithm described in Tarjan[Ref. 15], to reorder the ver- 
tices of the subgraph H into a new perfect elimination ordering. The computational 
complexity of this algorithm is O(n+m). The function mes.m uses a list x of size |V| 
representing each of the vertices, which 1s initially set to zeros. Each time a vertex 
visited all of its neighbors the value on list zx is increased by one. We use the vertex 
with the highest number on the list x as our next vertex to visit. We add this vertex 
to the list v and then continue the search. The algorithm terminates when all the 
vertices have been added to v. The list v represents the perfect elimination ordering 
of H. This new ordering ends with the first vertex of the old perfect elimination 


ordering, and is still a’perfect elimination ordering since the subgraph H 1s chordal. 


8. Edge Completion 

We now use a procedure first described by Grone, Johnson, et al in [Ref. 16] to 
perform a edge completion sequence on the chordal graph until all missing edges have 
been reinserted. This results in a supergraph of the input graph. The computational 
complexity of this procedure is O(n). In the function complete.m we start at the 
last vertex in the perfect elimination ordering and connect it to the remaining vertices 


starting with the next highest. Each time an edge in added to the graph we check to 


uy 


see if it is on the list of missing edges. If it is on the list of missing edges we eliminate 


it from the list and continue the procedure until the list is empty. 


9. Greedy Coloring 

Finally we perform a procedure on the chordal supergraph to determine its 
chromatic number with the function grcolor.m. This function takes advantage of the 
input chordal graph in a perfect elimination ordering. Since S is a perfectly orderable 
graph, applying the greedy coloring algorithm produces an optimal coloring of the 
graph in O(n). Each vertex is colored using the cheapest available color until all have 


been colored. 


C. IMPROVED ALGORITHMS 


With these improved algorithms we want to reduce the relative difference in 
the chromatic numbers of the original arbitrary graph and the chordal supergraph. 
First we experiment with manipulations on the arbitrary graph and then we examine 


some special ordering of the chordal graphs. 


1. Vertex Sort Algorithm 

This variation is based on the assumption that a better perfect elimination 
ordering is produced for the subgraph H when the vertices of G are pre-sorted by 
degree, highest first. This perfect elimination ordering might facilitate the selection of 
the supergraph S that would minimize the amount of added edges. In this variation we 
sort the vertices of the arbitrary graph G with the function versort.m. The vertices 
with the highest degree are considered first in the function mchord.m. Thus, the 
resulting maximal chordal subgraph H has an improved perfect elimination ordering. 


This algorithm is incorporated in the function project2.m. 


2. Missing Edge Algorithm 
This variation presupposes that a better perfect elimination ordering is used on 


the maximal chordal subgraph H if the vertices of H are presorted so that the vertices 


with the most missing edges are considered first in the edge completion algorithm. 


The resulting supergraph S should contain the minimal number of unnecessary edges. 


D. COMPUTATIONAL COMPLEXITY 


The most time-consuming function in this algorithm is the greedy-backtracking 
coloring function gbcolor.m, which is an inefficient non-polynomial time algorithm. 
All other functions in these algorithms run in polynomial time (see Table I). The 
inefficiency of gbcolor.m has limited us to graphs of order 20 or less. We presume 


another limiting factor in this algorithm is the programming language MATLAB. 


MATLAB Function | Worstcase Run-time 





unigraphi.m Olne) 

connecti.m O(n+™m) 

mchord.m O (mA) 

mclique.m O(n+m) 
gbcolor.m non-polynomial time 
mcs.m O(n +m) 
complete.m Cg) 

grcolor.m O(n) 


Table I. Computational Complexity. 


Zl 








IV. EXPERIMENTAL RESULTS 


The algorithms in chapter III were implemented in MATLAB and tested on 
100 random graphs of various densities (where density is the probability of an edge 
existing between any two vertices), having an order of 7 to 20. The projects were run 
on a HP700/15 workstation. Table II summarizes the results. The results are stated 
in the framework of the relative error between computing the actually chromatic 


number of the graph G and the chromatic number of the chordal supergraph S. 


A. GRAPH ORDER 
As can be seen from Table II and Figures 5 - 7, relative error increases with 
the order of the graph. Problem difficulty for our algorithms increases as well, which 


is not peculiar to our approach, but is intrinsic to the nature of the problem. Graphs 
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Figure 5. Relative error of y versus graph order. Note: This is discrete data. Relative 
error has been represented by continuous lines for clarity only. 
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Figure 6. Relative error of y versus graph order. Note: This is discrete data. Relative 
error has been represented by continuous lines for clarity only. 


of large order tend to be harder to color due to the increased complexity of the graph. 
We are limited to graphs of order 20 due to the excessive use of computer time. The 
computation time of the function gbcolor.m to compute the exact y(G)in graphs 
larger than order 20 was too excessive. The Missing Edge algorithm shows promise 
because the relative error is less than the other two algorithms even in the higher 
orders. On the other hand, the Vertex Sort algorithm does not perform as well in 
the lower orders as the other two algorithms, and does not show promise for further 
research. The Basic algorithm was used in this circumstance to provide a basis to the 


improved algorithms. 


B. GRAPH DENSITY 
The density of the graph is derived from the probability p that a given edge 
occurs. As can be seen from Table II and Figures 5 - 7, the density of the graph 
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Figure 7. Relative error of y versus graph order. Note: This is discrete data. Relative 
error has been represented by continuous lines for clarity only. 


has a peculiar effect on the performance of all the algorithms. The edge probability 
p = .2 produces the largest error in all three algorithms at the higher order, and this 
error appears to be increasing faster than when other probabilities are used. There 
appears to be a peculiarity in the relative error of y for graphs with edge probability 
between .1 and .3. Further research on the structure of random graphs with these edge 
probabilities is required to shed some light on this peculiarity. Again, the Vertex Sort 
algorithm shows little promise since the characteristics of Figure 6 shows no distinct 
pattern of the relative error in y for different probabilities p. The Missing Edge 


algorithm shows promise especially in the probabilities p = .1 and p =.5. 


C. CONCLUSION 
In conclusion, we see that we have moved the NP-complete problem from the 


coloring of an arbitrary graph to the choice of the correct maximal chordal subgraph 
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with the perfect elimination order which will produce the minimum chordal super- 
eraph of the arbitrary graph. If it were possible to choose this maximal chordal 


subgraph accurately the relative error of y would be. 
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Table II. Experimental Results of the Supergraph Heuristic. 
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V. FURTHER RESEARCH 


We have given a new approach to coloring an arbitrary graph with the use of 
a supergraph heuristic. With higher order graphs the relative error of the estimate 
is larger than anticipated, however the procedure for computing the estimate is a 
polynomial-time algorithm. The Edge Sort preprocessing algorithm shows promise 
of improving the relative error of the estimate and should be further studied for 
continued improvement. The main area for further research will be in the choice of 
the maximal chordal subgraph which will produce the minimum chordal supergraph. 
If this choice can be performed accurately our relative error would be next to nothing. 

It follows from the work of Grone, Johnson, et al. that if G’ is any chordal 
supergraph of G, and if His the chordal subgraph of G produced by Maxchord or 
some related algorithm, then there exist a chordal completion sequence containing 
both H and G’. Finding such a sequence would be ideal. Since the algorithm for 
generating such sequences is completely driven by perfect elimination orderings, our 
problem reduces to that of finding an optimal perfect elimination ordering for H. 

Further study is required in the analysis of the relative error of the estimate 
for the data in Table II and for larger graphs. To acquire data from larger graphs 
we must compile our programs into a more efficient programing language, possibly 
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APPENDIX A. PROGRAM FOR BASIC 
ALGORITHM 


7, Loren Eggen Project #1 revised 30 May 97 


p=.1; 4 probability of edge present 
fid=fopen(’Results2/projecti.out’,’a’); 

iemintt(fid,’Project #1\n’); 

fprintf(fid,’Edge probability in Arbitrary Graph is:%44.1f\n’,p); 
fprintf(fid,’ Matrix Ave. Chromatic Numbers Elapsed\n’); 


Boeintt(fid,’ Size Actual New Error Time\n’) ; 

telose(fid) ; 

mer 1i-/:20 order of arbitrary graph G 
M=[]; initialize storage matrix 
k=0; initialize k 


h 

h 

yf 

t=clock; i> start clock 

while k < 100; 4 generate k graphs 
v=L]; 4, .imitialize v 
while length(v) “= i 4 loop for undirected connected graph 

G=unigraph1(i,p); 4 generate a random undirected graph G 

[v,e]=connecti(G,1); 4, check if G is connected 

end; %, end while 

{v1,ei1]=mchord(G,1); 4 find the maximal chordal subgraph H 

H=mkadjmat(e1,i); 4 make adjancency matrix for H 

F=G-H; 4 determine missing edges 

em=mkedges(F); clear F; 4 make a list of missing edges 

if ~isempty(em) ; 4 if the original graph was not chordal 
Fo HGviled ie lear H- %4 reorder vertices of H to mchord peo-order 
emi=trans(vi,em) ; 4 translate missing edges to new order 
c1=mclique(H1) ; % maximum clique of Hi, lower bound of G 


vcl=gbcolor(G,c1); 4 color G using greedy-backtracking coloring 
ac=max(vci); clear G; 7, chromatic number for G 

v2=mcs(H1,1); 4 find peo ordering using max. card. search 
H2=H1i(v2,v2); clear H1: 4 reorder vertices Hi to mcs peo-order 
em2=trans(v2,em1) ; 7, translate missing edges to new order 
S=complete(H2,em2) ; 4 complete H2 until missing edges are added 
c2=mclique(S); clear H2; 4 maximum clique of S 
v=fliplr(1i:length(S)); 4 reverse the ordering of S 
vc2=grcolor(S(v,v) ,c2); 4 color S using greedy coloring 
nc=max(vc2); clear S; 74, chromatic number for S$ 

M=(M;i ac nce (nc-ac)/ac]; 4, record results 
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h=k+1 : 4, increment k 


end; %4 end if 
end; % end while 
[n,m]=size(M) ; 4 size of M 
aac=sum(M(:,2))/n; f, average actual chormatic number 
Sic=sumiMi- =o) n- 4, average new chormatic number 
adc=sum(M(:,4))/n; /,, average relative error 
Eim—ectime vclock.t)/ oe /, elapsed time 


fid=fopen(’Results2/project1.1.out’,’a’); 
fprintf(fid,’Edge probability in Arbitrary Graph is:/4.1f\n’,p); 
fp ide (tad eiMantiais Chromatic Numbers\n’) ; 
fpRintt tidy waarec Actual New Error\n’); 
Pomintt £id, ~orOtmnn Of S/oimO tern a eee): 
fprintf(fid, ’Average:/5.0f 47.0f 47.4f\n’ ,aac,anc,adc); 
felloce Guid: 
fid=fopen(’Results2/projecti.out’,’a’); 
fprintt (fid,°/5.0f (7e0r 47. 0f 47 .4f (8. 1f\n’ ,1,aac,anc,adc, time 
helocet fide 
end; 4 end for 
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APPENDIX B. PROGRAM FOR THE 
VERTEX SORT ALGORITHM 


7, Loren Eggen Project #2 


=a 


revised 30 May 97 


4 probability of edge present 


fid=fopen(’Results2/project2.out’,’a’); 


apeintt(fid,’ Project #2\n’); 


fprintf(fid,’Edge probability in Arbitrary Graph is:%4.1f\n’,p); 
fprintf(fid,’ Matrix Ave. Chromatic Numbers Elapsed\n’); 


memintt(fid,’ Size Actual New 


fclose(fid) ; 
for 1=7:20 
M=(] ; 
k=0; 
t=clock; 
while k < 100 
v=([]; 
while length(v) “=i 
G=unigraphi(i,p); 
[y,el=connect1(G,1); 
end; 
[G1 ,vi=versort(G): 
[vi,ei]=mchord(Gi,1); 
H=mkadjmat(ei,i); 
F=G1-H; clear G1; 
em=mkedges(F); clear F; 
if ~isempty(em) ; 
Hi=H(vi,vi); clear H; 
emi=trans(vi,em); 
c1=mclique(H1) ; 
vct=gbcolor(G,c1) ; 
ac=max(vci); clear G; 
v2=mcs(H1,1); 
H2=H1i(v2,v2); clear H1; 
em2=trans(v2,emi); 
S=complete(H2,em2) ; 
c2=mclique(S); clear H2; 
v3=fliplr(i:length(S)) ; 
ViG2moueolor(S(v3.,Va),c2) ; 
nc=max(vc2); clear S; 
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Error Time\n’); 


order of arbitrary graph G 

initialize storage matrix 

initialize k 

start clock 

number of graphs to generate 

initialize v 

loop for undirected connected graph 
generate a random undirected graph G 
check if G 1s connected 

end while 

sort the vertices by highest degree 

find the maximal chordal subgraph H 

make adjancency matrix for H 

determine missing edges 

make a list of missing edges 

if the original graph was not chordal 
reorder vertices of H to mchord peo-order 
translate missing edges to new order 
maximum clique of Hi, lower bound of G 
color G using greedy-backtracking coloring 
chromatic number for G 

find peo ordering using max. card. search 
reorder vertices Hi to mcs peo-order 
translate missing edges to new order 
complete H2 until missing edges are added 
maximum clique of S 

reverse the ordering of S 

color S using greedy coloring 

chromatic number for $ 


Melitta ae Wemcnc—ac) /acie 4, record results 


k=k+1; 4% increment k 
end; 4, nd if 

end; fmendewhabie 
[n,m]=size(M); 4, average actual chormatic number 
aac=sum(M(: ,2))/n; 4 average new chormatic number 
anc=sum(M(:,3))/n; ,, average relative error 
adc=sum(M(: ,4))/n; 4 average relative error 
tim=etime(clock,t)/60; 7 elapsed time 


fid=fopen(’Results2/project2.1.out’,’a’); 
fprintf(fid,’Edge probability in Arbitrary Graph is:%4.1f\n’,p); 
foment n hid. enianta x Chromatic Numbers\n’) ; 
Epuintcl tid,’ sysaze Actual New  Error\n’); 
Epigint hid.) ,SfO0LageOl AO 1,7 <4e me ,M?) ; 
fprintf(fid,’Average:/5.0f 47.0f “47.4f\n’ ,aac,anc, adc); 
helose (tad) 
fid=fopen(’Results2/project2.out’,’a’); 
Fomine (tid, yo.0f 47.02 4.0L 47 .4f 48.1f\n’ ,1,aac,anc,adeyuumme 
Helose (rfid) 
end; 4 end for 
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APPENDIX C. PROGRAM FOR MISSING 
EDGE ALGORITHM 


/, Loren Eggen Project #3 


p=.1; 


revised 30 May 97 


4, probability of edge present 


fid=fopen(’Results2/project3.out’,’a’); 


momint?(fid,’*Project #3\n’); 


fprintf(fid,’Edge probability in Arbitrary Graph is:44.1f\n’,p); 
fprintf(fid,’ Matrix Ave. Chromatic Numbers Elapsed\n’); 


Pprinue(fid,? Size Actual New 


fhelose(fid); 
for 1=7:20 
M=[] ; 
k=0; 
t=clock; 
while k < 100 
v=(]; 
while length(v) “=i 
G=unigraph1(i,p); 
[v,e]=connect1(G,1); 
end; 
[v1,el1J=mchord(G,1); 
H=mkadjmat(e1,i); 
F=G-H; 
em=mkedges(F); clear F; 
if ~isempty(em) ; 
H1i=H(v1,vi); clear H; 
emi=trans(vi,em); 
c1=mclique(H1) ; 
vcil=gbcolor(G,c1); 
ac=max(vc1); 
G1=G(vi,vi); clear G; 
F1=Gi-H1; clear Gi; 
Ly, vs]=sort(sum(F1)); 
v2=fliplr(vs); clear F1; 
H2=Hi(v2,v2); clear H1; 
em2=trans(v2,em1); 
v3=mcs(H2,1); 
H3=H2(v3,v3); clear H2; 
em3=trans(v3,em2) ; 
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Error Time\n’); 

order of arbitrary graph G 
initialize storage matrix 
initialize k 

start clock 

number of graphs to generate 
initialize v 

loop for undirected connected graph 
generate a random undirected graph G 


check if G is connected 

end while 

find the maximal chordal subgraph H 

make adjancency matrix for H 

determine missing edges 

make a list of missing edges 

if the original graph was not chordal 
reorder vertices of H to mchord peo-order 
translate missing edges to new order 
maximum clique of H1, lower bound of G 
color G using greedy-backtracking coloring 
chromatic number for G 

translate G into H1 ordering 

find the max. edges missing 

sort by max. edges missing 

descending order 

sort Hl max. edges missing first 
translate missing edges to new order 

find peo ordering using max. card. search 
reorder vertices H1 to mcs peo-order 
translate missing edges to new order 


S=complete(H3,em3) ; 4 complete H2 until missing edges are adde 


c2=mclique(S); clear H3; 4, maximum clique of § 
v4=fliplr(1:length(S)) ; 74 reverse the ordering of S 
vc2=grcolor(S(v4, v4) ,c2) ; 4 color S using greedy coloring 
nc=max(vc2); clear S; 4, chromatic number for S$ 
M=[(M;i ac nc (nc-ac)/ac]; 4, record results 
k=k+1; /, increment k 
end; 7, endear 
end: 4, end while 
[n,mJ=size(M); 4, size of M 
aac=sum(M(: ,2))/n; /, average actual chormatic number 
anc=sum(M(: ,3))/n; 4 average new chormatic number 
adc=sum(M(: ,4))/n; /, average relative error 
tim-ebtime clock, t)/60- 4, elapsed time 


fid=fopen(’Results2/project3.1.out’,’a’); 
fprintf(fid,’Edge probability in Arbitrary Graph is:%4.1f\n’,p); 
seyene allo neae (Caealol wensrest ¥< Chromatic Numbers\n’) ; 
forint merad, Size Actual New  Error\n’); 
ipigimiciirad, 7,5. Of oi nOL 4.7 OLe),/ .4f\n’? M7) ; 
Meraner (tid, Average:,5,.01 4/.0f 47.4f\n’ ,aac,anc, adc) ; 
Polose (Lid); 
fid=fopen(’Results2/project3.out’,’a’); 
PoGeiInvec tdeewo. Of (Oke nor ,/ of *sS 1f\n’,1,aac,ane, ade .campe 
felose (fu): 

end; 4, end for 
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APPENDIX D. GRAPH COMPLETION 
FUNCTION 


function A=complete(A,e1) 


4 function A=complete(A,e1) 
y/ 
4 This function is a graph completion function for the thesis project. 
4 Input a peo ordering v, a list of edges e from a maximal chordal subgraph, 
i, and a list of edges e1 necessary to make a super-hypergraph of the original 
4 graph. The output is a chordal supergraph. 
4 by Loren G. Eggen, 14 April, 1997. 
A=A+speye(size(A)); 7, add loops 
[n,m]=size(A); /eteot veruices 
vei:in; /pliseeOLeventices 
p=sum(A); q=find(p~=n) ; 74 eliminate and full 
v=(v(q)); 4 vertices from list v 
while ~isempty(e1) 4 loop for all missing edges 
k=length(v) ; 4 vertix by peo ordering 
l=max(find(~A(v(k),:))); 4 next highest missing edge 
if ~isempty(1) 4 test if found next 
ez=(v(k) 1]; /, edge to add ez 
e=[e; ez]; i, add edge to list 
A(v(k) ,1)=1; / add edge in adjacency 
Alek) )=1- emote 
end; % end of if 
if sum(A(v(k),:))==n % test if vertix full 
h=find(v~=v(k)); 7, eliminate vertix from 
v=v(h); 7, Ist of vertices 
end; i endeon at 
[a,b]=size(e1) ; 4 a = length of missing edges 
for i=l:a 74 for each missing edge 
if all(e1i(i,:)==ez | e1(i,:)==fliplr(ez)) 4, test ez in missing edges 
if a== 4 if yes and last then 
e1=(]; 4 empty list of missing edges 
break; 4 break for loop 
else i Deonobk last 
x=1:a; % set x 
x=x(find(x7=i)); 7, find missing edge = ez 


oT 


el=el(x’,:); 7, delete edge from list 


break; 4, break for loop 
end; ) cudennered + 
end; Peondeol cer sit 
end; 4 end for loop 
end; end ot whale 
A=A-speye(size(A)); 4 eliminate loops 
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APPENDIX E. CONNECTED GRAPH 
FUNCTION 


function [v,e]l=connecti(A,i); 


4 function [v,e]=connect1(A,1i) ; 


h 


4 This function will find a connected component in‘the input graph A and 
4, starting vertex 1. It uses depth first search and outputs the vertices 
4 and edges of the connected component. 


74, By Loren G. Eggen, 30 May, 1997. 


V=l; 
v0=1; 
e=(]; 
n=length(A) ; 
while length(v) <n 
x=adj(A,1); 
if ~isempty(x) 
t=0; 
for j=1:length(x) 
if isempty(find(x(j)==v)) 
velv,x(Qj)J; 
e=[e; i,x(j)J; 
BCs eG) )=0; 
A(x(j) ,1)=0; 
i=x(j); 
t=1; 
break; 
end; 
end; 
os Gea 
for j=1:length(x) 
e=[e;i,x(j)J; 
ACi,x(j))=0; 
ACO). =0; 
end; 
l=1; 
while sum(A(i,:)) == 0 & 


initialize list of vertices 
first vertex 

initialize list of deleted edges 
number of vertices 

loop till all vertices are added 
adj. vertices to current vertex 
if x is not empty 

test variable if vertex 1s added 
for all the adj. vertices 

find 1st one not inelist Vv 

add it to the list v 

update deleted edges 

eliminate edge in adj. matrix 
eliminate edge in adj. matrix 
make new vertex current 

set test variable true 

break when new one found 

end of if 

end of for 

if no new vertex but x not empty 
for each adj. vertex which is on v 
update deleted edges 

eliminate edge in adj. matrix 
eliminate edge in adj. matrix 
end of for 

set backtracking index 
backtrack till edge is present 


1=v(length(v)-1); 
L=i+1 ; 
end; 
end; 
elseif 1 == v0 
break; 
else 
1l=1; 
while sum(A(i,:)) == 0 & i ~= v0 
i=v(length(v)-1); 
l=l1+1; 
end; 
end; 
end; 
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backtrack list v 

increment index 

end of while 

end of if 

1f x was empty and we returned to v0 
break while loop, output component 
not at the start but x is empty 
set backtracking index 

backtrack till edge is present 
backtrack list v 

increment index 

end of while 

end of if 

end of while 


APPENDIX F. GREEDY-BACKTRACKING 
COLORING FUNCTION 


function Lvc]=gbcolor(A,m) 
4, function [vc]=gbcolor(A,m) 


4 This function uses the greedy-backtracking approach to 

4 color the vertices of a graph so that no two colors are 

4, together. Input the adjacency matrix of the graph and a 

4 minimum number of colors. Output ve is the vector of vertix 
colors. 


4 by Loren G. Eggen, 18 March, 1997. 


i=1; 4 starting index 
v=i; 7, first vertex 
vc=i1; /oGitesteu en bens scolor 
n=length(A); 4 number of vertices 
while length(vc) < n; 4 used until all vertices have been colored 
k=1; J Grrst color 
i=i+1; 4, increment index 
v=(v i]; JevectomiGh vertices Visited 
x=adj(A(v,v),i); 4 adjacent visited vertices of index 
xc=sort (vc(x)); 7, sorted colors 
for j=1:length(xc) 4 find the next available color 
if xc(j) == 4 if current color used 
k=k+1; yeinecrement color 
end; 4 ema if 
end; 7%, end for 
7 Kk > m Apt We ain out of Golors backtrack 
i=i-1; 7, decrement index 
vev(1:1); 4 go back one vertix 
ve=ve(1:i); 4 eliminate last color if necessary 
while length(v) > 1 4 do not backtrack past 1 
t=0; 4, test variable to break backtrack 
while vce(i) <m 4, if the color is < max see if we can 
VeECip=yoq)ti- 4, increase the color 
if “any(vc(i) == vc(adj(A(v,v),i))) 4, test if the color has been 
t=1; break; 4 used if not use it and stop 
end; yA) Xe Mb a 


4] 


ve=ve l:1); 

end; 

end; 

if length(vc) < length(v) 
ve=([vc kl]; 

end; 

if length(v) == 1&t == 
m=mt+i; 

end; 

end; 
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end while k < m 

found one that could be incremented 
break outer loop 

end if 

decrement index 

go back one vertix 

eliminate last color if necessary 
end while v > l 

end if k > m 

ifecnewinner loop did not 

update the color 

end if 

if we have ran out of colors and 
backtracked to the origin increase 
the available colors 

end while vc <n 


APPENDIX G. GREEDY COLORING 
FUNCTION 


function [vc]=grcolor(A,m) 


4 function [vc]=grcolor(A,m) 


4 This function uses the greedy approach to color the vertices of a graph 


4 so that no two colors are together. 


Input the adjacency matrix of the graph 


% and a maximum number of colors. Output vc is the vector of vertix colors. 
4 Optimal coloring if the input graph is chordal and reverse order perfect 


% elimination scheme. 
4 by Loren G. Eggen, 23 April, 1997. 


n=length(A) ; 


k=1; 
1=1; 
V=l1; 
vc=k; 
while length(v) < n; h 
k=1; i 
i=it1; ts 
v=[v il; i 
x=adj (A(v,v),1i); h 
xc=sort (vc(x)); ib 
for j=1:length(xc) h 
if xc(j) == y/ 
k=k+1; h 
end; i 
end; if 
vo=[vc k]; if 
if max(vc) > m if 


fprintf(’colors used greater than 
end; 
end; 
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used until all vertices have been visited 
mIrst color 

increment index 

vector of vertices visited 
adjacent visited vertices of index 
sorted colors 

find the next available color 

if color used 

incremet color 

end if 

end for 

update the color 

if we have ran out of colors 
Colors eiven wun); 

end if 

end while 





(f 
oon a 
x 7 
a _ 
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APPENDIX H. MAXIMAL CHORDAL 
SUBGRAPH FUNCTION 


function [v, E] = mchord(A,i) 
4 function [{v, E] = mchord(A,i) 


4, Returns the peo ordering of vertices and a set of edges 

4 which will generate a maximal chordal subgraph. Adjacency 
4 matrix A should represent a connected undirected graph. 

% This function uses algorithm MAXCHORD, P.M. Dearing. 


4 by Loren G. Eggen, 6 February, 1997 
% adj.m function called 


4 begin mchord 


v=1; 4, starting vertex v 
n=length (A) ; 74 number of vertices 
=i ie Vmibrsteotavertices 
s=s(find(s~=i)); 74, delete first vertex from the list 
y/ 
yf 
if 


C=zeros(n) ; initialize the set of adj. vertices 
C=sparse(C) ; 
E=(] ; initialize set of edges 
fom j—1:n-1 4 loop through all vertices except ist 
[vi e1]=adj(A,i); 4 find adj. vertices to current 
for k=1:length(v1) 4, loop through each adj. vertex u 
mo stmamecuiik) 7: )"l=C(i.,:)* ) face sctmGanmcibser sor C(y) 


make matrix sparse 


iietest.==—-C(i,-) 4 if so then 
Civ i(k), 1) =1; 4 update C(u) 
F=(E;i vitk)]; 74, update set of edges 
ACi,viCk))=0; 4 delete edge from adjacency matrix 
A(vi(k) ,i)=0; 7, both edges 
end; jh avol Ope ste 
end; / end@ote inne: for 
[1 m]=max(sum(C(s,:)’)); 4 find next vertex with largest card. in C 
v=[v s(m)]; 4 assign new vertex to reverse peo ordering 
i=s(m) ; 4, assign new vertex to current vertex v 
s=s(find(s~=i)); Jadebete vetlom list of aventices 
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end; Zoend Of ollcer for 
v=flipir(v) ; 4, peo ordering 


4, end of mchord 
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APPENDIX I. MAXIMUM CLIQUE 
FUNCTION 


function c=mclique(H) 


% function c=mclique(H) 


4% This function calculates the maximum clique number of a triangulated 
4 graph H, which is ordered by it perfect elimination scheme. 


74 by Loren G. Eggen, 21 April, 1997 


4 calls adj.m function 


c=1; 
{n,m]=size(H); 
S=zeros(1i,n); 
a=1:n; 
for 1=1:n 
X=(]; 
v=a(i); 
[x,e]=adj(H,v); 
for j=1:length(x) 
if find(v==a) < find(x(j)==a); 
Kx): 
end; 
end; 
if isempty(x), v; end; 
if ~isempty(X) 
u=min(X) ; 
S(u)=max(S(u) , length(X)-1) ; 
if S(v) < length(X) 
lei 3] 
c=max(c,1+length(X)); 
end; 
end; 
end; 


A’ 


initial clique number 

order of input graph 

list, # if times vertices i:n visited 
vertices 1:n 

loop for each vertex 

initialize X 

assign v current vertec 

find adj. vertices to v 

loop for each ad}. vertex 

if index of v < index of adj. vertices add 
adj. vertex to the list of higher indices 
end if 

end innner for 

if no adj. vertices v is cluque 

if X not empty 

u, smallest index in X 

assign S(u) max. of current value or clique 
if S(v) < current clique 

print current clique 

update maximum clique number 

end inner if 

end outer if 

end for loop 
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APPENDIX J. MAXIMUM CARDINALITY 
SEARCH 


function [v] = mcs(A,i) 

¥ function [v] = mcs(A,i) 

4 Returns a vector of vertices which indicate 

4 a possible perfect elimination scheme. Adjacency 

%Z matrix A should represent a connected undirected 

4 graph. This function uses Maximum Cardinality Search. 
4 by Loren G. Eggen, revised 29 January, 1997 


4 adj.m function called 


2, begin mcs 


v=i; 4 1st vertex in the peo 

n=length(A) ; 4 number of vertices 

x=ones(1,n); 7% initialize cardinality vector x 

for j=i:n-1 4, loop for each vertex 
vizadj(A,i); 7, find adj. vertices vi to current vertex 
x(v1)=x(v1)+1; 4 update the cardinality of vertices in vi 
[k,1]=max(x) ;: 4 find the vertex in x with max. cardinality 
velv 1]; 7, add new vertex to the peo list v 
x(v)=x(v)-x(v); i2Zero thesentricssoL xsster vertices inv 
A(i,1)=0; 7, eliminate edges from adj. matrix 
A(1,1i)=0; 74 eliminate edges from adj. matrix 
1=1; 4 assign current vertex to the last vertex added 

end; 4 end for loop 

v=fliplr(v) ; 7, reverse peo ordering 


4, end of mcs 
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APPENDIX K. MISCELLANEOUS 
FUNCTIONS 


function [v, e] = adj(A,i) 

4% function [v, e] = adj(A,i) 

4 Returns a vector of adjacent vertices and a list of 
i edges to the vertex i from the adjacency matrix A. 
4 by Loren G. Eggen, 29 January, 1997 

4% no intrinsic functions called 


4 begin adj 


v=]; (eanit Veieize v 
see i sbochesle wl 7 Aaee 
for j=1:length(A) 4, for each element in row i 
een 1, j) if an edge exist 
velv jl; 4 update v, list of adjancect vertices 
e=[e;i jl; 4 update e, list of associated edges 
end; / eenoent f 
end; 4, end for 
4 end of adj 


function A=mkadjmat(e,v) 

4 function A=mkadjmat(e,v) 

4 Input a set of edges containing the numeric vertices e.g. 
4 e=[1 2;1 7;2 3;2 5;3 4;4 6;5 6;6 7] and maximum vertices v. 
% Output an adjacency matrix of unidirectional graph. 


4 by Loren G. Eggen, revised 3 May, 1997 


A=zeros(v,v); 
([n,m]=size(e); 


O1 


£OG 11m 
ACeCig 1) O01, 27a 
Ge nee Gingell, 

end ; 

A=sparse(A) ; 

function e=mkedges(A) 

4 function e=mkedges(A) 


% This function makes a set of edges from the given 
7% adjancency matrix. 


4 by Loren G. Eggen, 11 April, 97 
K—-triwch). 
for i=1:length(A) 
y=rindaeaci-)) 
e=[e;i*xones(length(y) ,1),y]; 
end; 
function el=trans(v,e) 
4, function el=trans(v,e) 
,/ This function translates the edges into the new vertix ordering. 


4 by Loren G. Eggen, 14 April, 1997. 


[n,m]=size(e); 
Ory i=in 


ei(i,1i)=find(e(i,1)==v); 4 translate the missing edges 
e1(i,2)=find(e(i,2)==v) ; 4 into the new ordering 
end; 


function [A] = unigrap1(n,p) 
4 function[A] = unigrap1i(n,p) 


4 Generates an edge with probability p in an adjacency 
A matrix for a unidirected graph. 
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% by Loren G. Eggen, revised 29 May, 1997 


”% no intrinsic functions called 


4 begin 

A=rand(n) ; 4 generate random 0-1 matrix nxn 
A=A(:,:)<p; 4 eliminate all entries > p 
A=A-diag(diag(A)) ; 4 eliminate diaganol 

A=triu(A); 4 eliminate lower triangular 
A=A+A’; 4 make symetric 

A=sparse(A) ; 4 make matrix sparse storage 

/% end 


function [S,v]J=versort (A) ; 
% function [S,v]=versort(A); 


4 Label an adjacency matrix sorting by highest degree 
feevertix in the matrix. 


i=sum(A) ; 

Ly, jJ=sort(i); 
v=fliplr(j); 
S=A(v,v); 


D3 
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